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1  Final  Technical  Progress 

1.1  Research  Conducted  by  Michael  Meyers 

Most  of  the  research  has  been  conducted  by  Mike  Meyers  developing  the  new  micro¬ 
physical  scheme  which  was  implemented  into  RAMS.  This  new  scheme  allows  prediction 
on  both  the  concentration  and  mixing  ratio  of  the  distribution  function  for  each  hydrom¬ 
eteor  category  making  the  determination  of  each  hydrometeor  size  spectra  less  arbitrary. 
With  more  detailed  physics  the  model  should  have  the  potential  to  allow  for  the  evolution 
of  the  hydrometeor  spectra.  The  new  microphysics  also  allows  prediction  of  an  additional 
hydrometeor  species,  hail.  Some  highlights  of  the  model  include: 

•  the  use  of  a  generalized  gamma  size-spectrum  where  the  v  parameter  can  be  prescribed 
by  the  user  as  opposed  to  a  fixed  Marshall  Palmer  spectrum 

•  the  introduction  of  ice-liquid  mixed  phase  graupd  and  hail  categories  categories  with 
non-thermal  equilibrium  for  the  rain,  graupd  and  hail  dasses 

•  new  heterogeneous  and  homogeneous  ice  nudeation  parameterizations 

•  approximate  solutions  to  the  stochastic  collection  equation  rather  than  the  continuous 
accretion  modd  approach 

•  breakup  of  rain  droplets  is  formtilated  into  the  collection  effidency 

•  analytical  flux  equations  predict  mixing  ratio  and  number  concentration  conversion 
from  pristine  ice  crystals  to  snow  due  to  deposition  (no  riming) 

•  predictive  equations  for  ice  nudd  (IN) 

•  crystal  habit  is  diagnosed  dependent  on  temperature  and  saturation 

•  evaporation  and  mdting  of  each  spedes  assumes  that  the  smallest  partides  completdy 
disappear  first 

•  more  complex  shedding  formulations  which  take  into  account  the  amount  of  water 
maas  on  the  coalesced  hydrometeor. 

After  testing  the  new  scheme,  two  case  studies  were  examined,  a  wintertime  orographic 
predpitation  event  over  the  Front  Range  of  Colorado  from  the  Winter  Idng  and  Storms 
Project  (WISP91),  and  a  summertime  convective  case  over  the  High  Plains  of  Montana  from 
the  Cooperative  Convective  Predpitation  Experiment  (CC0PE81).  These  cases  allowed 
comparisons  between  the  two-moment  predictive  scheme  and  the  one-moment  predictive 
scheme  and  highlighted  the  advantages  of  independently  predicting  on  two  moments  of 
the  hydrometeor  spectra.  Detailed  analysis  was  also  done  on  these  two  cases  with  spedal 
attention  given  to  radar  reflectivity,  aircraft,  and  ground-based  measurements. 

Other  research  conducted  during  this  period  indudes: 
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1.  quantitative  precipitation  forecasts  of  wintertime  orograpMc  cloud  from  the  Sierra 
Nevada  range  of  California 

2.  implementation  of  natural  and  artificial  primary  ice  nudeation  parameterizations  for 
use  in  doud  models 

3.  formulation  of  parameterizations  for  homogeneous  freezing  of  doud  droplets  and  haze 
partides  which  are  relevant  to  ice  initiation  of  cirrus  douds. 

1.2  Research  Conducted  by  Scot  C.  R.  Rafkin 

Mr.  Paflfin  (formerly  Randell)  has  been  funded  under  the  current  AASERT  proposal  for 
approximately  the  last  six  months.  His  research  has  been  directed  towards  the  development 
of  a  convective  parameterization  valid  for  numerical  models  with  grid  spacings  ran^ng  from 
the  mesoscale  to  GCM-like  scales.  A  summary  of  his  research  is  presented  bdow. 

1.2.1  Objective 

Currently,  a  cumulus  parameterization  which  is  valid  at  the  mesoscale  does  not  exist. 
A  few  schemes  have  been  developed  which  can  be  applied  at  scales  larger  than  about 
50-100km  {e.g.  Eritsch  and  Chappd  (1980);  Krdtzberg  and  Perkey  (1976)).  However, 
more  often  than  not,  parameterizations  whidi  have  been  devdoped  for  large-scale  modds, 
such  as  Arakawa  and  Schubert  (1974)  and  Kuo  (1965),  are  applied  at  smaller  mesoscales 
despite  the  approximations  and  assumptions  which  become  increasingly  questionable  as  the 
grid  spacing  decreases.  Furthermore,  modellers  are  faced  with  havmg  to  choose,  somewhat 
arbitrarily,  different  parameterizations  for  different  scales.  The  goal  of  this  research  is  aimed 
at  devdoping  a  cumulus  parameterization  which  is  vahd  over  a  wide  range  of  scales,  with 
particular  focus  on  at  the  mesoscale.  The  new  parameterization  will  eliminate  the  need  for 
different  parameterizations  at  different  scales,  and  free  the  modeller  from  having  to  dedde 
when  to  change  from  one  parameterization  to  another. 

1.2.2  Philosophical  Foundation 

Before  the  parameterization  can  be  designed,  we  must  consider  what  it  means  to  param¬ 
eterize  cumulus  convection  at  various  scales.  In  the  case  of  large-scale  parameterization, 
the  statistical  effects  of  an  ensemble  of  convection  is  considered.  We  define  an  ensemble  of 
douds  to  be  the  set  of  all  possible  douds  which  could  exist  in  a  given  environment.  In  the 
Arakawa-Schubert  parameterization  for  example,  all  possible  douds  which  could  exist  (the 
ensemble),  do  exist;  that  is,  they  contribute  to  the  cumulus  fordng  of  the  large-scale  envi¬ 
ronment.  As  long  as  the  area  we  are  considering  is  large  enough  to  contain  the  ensemble, 
this  approach  is  valid. 

Clearly,  as  a  modd  grid  box  shrinks,  a  point  will  be  reached  where  it  is  not  possible 
to  physically  fit  the  ensemble  within  the  given  area.  If  we  wish  to  know  the  effects  of 
cumlus  convection  on  the  modd-scale  variables,  however,  we  still  need  to  know  which  types 
of  douds  are  in  the  box.  We  define  the  active  population  to  be  the  subset  of  the  ensemble 
which  is  actually  found  within  a  ^ven  area.  We  can  make  an  educated  guess  at  the  active 
population  by  considering  which  type  of  douds  are  more  likely  to  be  found  within  a  given 
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area.  Of  course,  the  actual  cloud  population  within  the  area  will  generally  differ  from  our 
guess,  but  on  the  average  we  will  be  correct. 

Determining  a  probability  density  function  for  an  ensemble  of  clouds  will  be  one  com¬ 
ponent  of  the  parameterization.  Once  we  know  the  probability  for  the  existence  of  each 
type  of  doud  in  a  given  environment,  we  can  select  the  active  population  for  a  given  area. 
Large-scale  models  with  large  grid  spacings  will  contain  the  entire  ensemble  since  there  will 
be  room  for  even  the  most  unlikely  members.  As  the  grid  spacing  decreases,  we  will  use  the 
PDF  to  select  the  active  population. 

In  addition  to  the  selection  of  an  active  population,  we  must  also  consider  that  the 
effects  of  the  convection,  in  particular  the  subsident  circulations,  may  extend  beyond  the 
grid  box  boundaries.  In  large-scale  parameterizations  the  effects  of  the  cumulus  ensemble 
are  assumed  to  remain  wholly  within  a  grid  box — a  horizontally  local  event.  Therefore,  as 
we  transition  from  large  to  small  grid  spacings,  we  also  transition  from  a  horizontally  local 
parameterization  to  a  horizontally  non-local  parameterization.  We  conclude  that  the  area 
over  which  the  parameterization  is  applied  should  not  be  a  function  of  the  grid  resolution. 
Rather,  it  should  be  applied  over  the  convective  influence  area  (CIA),  or  the  area  over  which 
the  convective  circulations  extend.  Determining  the  CIA  will  be  another  component  of  the 
parameterization. 

1.2.3  Parameterization  Design 
Prognostic  Cumulus  Kinetic  Energy 

Following  the  work  of  Lord  and  Arakawa  (1980),  and  Randall  and  Pan  (1993),  we 
introduce  the  prognostic  cumulus  kinetic  energy  (CKE)  equation: 

^^  =  (A(A)-D(A))Mb(A),  (1) 

where  K  and  Mb  are  the  vertically-integrated  kinetic  energy  and  the  cloud-base  mass  flux 
per  unit  area,  respectively,  associated  with  the  sub-ensemble  of  clouds  in  the  interval  (A,  A+ 
dX).  The  doud  work  function.  A,  and  D  are  the  generation  and  dissipation,  respectively,  of 
K  per  unit  doud-base  mass  flux.  As  proposed  by  Xu  (1991)  and  Arakawa  and  Xu  (1990) 
we  rdate  K  and  Mb  through  the  constant  a  as  follows: 

K  =  aMl.  (2) 

Combining  Eqs.  (1)  and  (2),  and  assuming  that  a  is  independent  of  time  yields  a  linear 
differential  equation  for  the  time  rate  of  change  of  the  sub-ensemble  doud-base  mass  flux: 

OMb  _ 

dt  ~  2a  2Td  ’ 

where  the  dissipation  term  has  been  modded  as  D  =  K/vd-  Given  the  parameters  a  and  Td 
for  each  sub-ensemble,  we  can  predict  the  evolution  of  their  assodated  doud-base  mass  flux. 
Currently  we  are  using  fixed  values  for  both  parameters,  but  future  research  is  intended  to 
provide  more  realistic  parameterizations  of  these  quantities. 

The  above  method  of  determining  the  doud-base  mass  flux  is  in  contrast  to  the  dassical 
technique  typically  employed  in  the  Arakawa-Schubert  parameterization.  We  note  that  the 
prognostic  equation  pven  by  Eq.  (3)  is  computationally  much  more  effident.  In  addition,  it 
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does  not  require  partitioning  or  separating  the  large-scale  forcing  from  the  cumulus  forcing. 
The  implications  of  the  closure  are  discussed  in  more  detail  in  Randall  and  Pan  (1993). 
Determitiing  the  PDF  and  Active  Population 

The  cloud  work  function  for  each  sub-ensemble  is  computed  using  a  plume  model.  Fol¬ 
lowing  Simpson  and  Wiggert  (1969)  we  can  calculate  the  in-doud  vertical  vdodty  proiile 
and  kinetic  energy.  Let  K(X)  be  the  vertically  integrated  kinetic  energy  of  the  plume-doud 
of  t3rpe  A.  Then, 

KiX)A  =  n(A)o(A)2i:,  (4) 

where  A  is  the  area  in  which  the  sub-ensemble  is  present,  n  is  the  number  of  A  douds  and 
a  is  the  vertically  averaged  cross-sectional  area  of  the  a  type  A  doud.  The  quantity  nc 
is  the  total  area  occupied  by  each  sub-ensemble,  and  is  used  to  determine  the  probabiUty 
distribution  function  (PDF)  and  active  population.  By  normalizing  this  quantity,  we  obtain 
the  PDF. 

There  are  two  methods  to  determine  the  active  population.  The  first  method  involves 
rlinnBiTig  the  most  probable  doud  types  from  the  PDF,  and  allowing  only  as  many  sub¬ 
ensembles  that  will  fit  within  the  modd  grid  box.  The  total  area  of  the  active  population 
is  simply  the  sum  of  no  over  the  chosen  doud  types.  This  procedure  will  only  provide 
the  effects  of  the  most  probable  douds  and  never  indude  any  effects  from  more  unlikdy 
doud  types.  After  choosing  the  active  population,  we  set  the  doud-base  mass  flux  for  each 
sub-ensemble  not  sdected  equal  to  zero.  Of  course,  this  allows  the  doud  work  function  for 
the  unsdected  sub-ensembles  to  increase  with  time  since  there  is  no  stabilizing  convective 
effects  allowed.  Eventually,  the  cloud  work  will  increase  enough  that  the  douds  “move  up” 
the  probability  scale  and  are  sdected. 

In  the  second  method,  the  quantity  na  is  normalized  by  the  grid  box  area.  This  new 
distribution  function  then  represents  the  fraction  of  each  sub-ensemble  that  could  exist 
in  a  grid  box.  The  doud-base  mass  fluxes  are  then  multiplied  by  this  quantity  when  de¬ 
termining  the  doud  fordng.  In  this  method,  all  douds,  even  the  unlikdy  ones,  have  an 
influence.  Following  the  determination  of  the  active  population,  the  doud-base  mass  fluxes 
axe  normalized  by  the  PDF. 

It  is  not  dear  yet  which  method  should  be  chosen.  Once  preliminary  testing  is  completed, 
we  win  undert2tke  simulation  using  both  methods  to  try  and  determine  the  more  appropriate 
method. 

Computing  Grid-scale  Tendendes  Due  to  Clouds 

Given  the  a<:tive  population  within  a  grid  box,  we  can  determine  the  average  tendendes 
of  the  environment  variables,  ^e>  d'l®  fo  each  sub-ensemble  in  the  active  population  through 
the  divergence  of  the  cumulus  fluxes  which  are  modeled  as: 

=  M(^c-  6)>  (5) 

where  M  is  the  total  doud  mass  flux  due  to  the  active  population,  and  the  subscripts  c  and 
e  refer  the  doud  and  environment  properties,  respectivdy.  However,  the  average  tendency 
is  not  valid  over  the  modd  grid.  It  is  valid  over  the  CIA  for  the  particular  sub-ensemble 
which  may  in  general  cover  many  modd  grid  boxes.  We  therefore  need  to  determme  the 
CIA  assodated  with  eadi  particular  sub-ensemble  within  the  active  population. 

We  assume  that  the  subsident  circulations  due  to  a  doud  are  bounded  by  the  wavefront 
of  an  internally  propagating  gravity  wave.  This  assumption  has  been  verified  by  observation 
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(Mapes,  1993)  and  also  by  numerical  experiments  (NichoUs  et  al,  1991).  The  speed  of  the 
wavefront  is  obtained  through  linear  gravity  wave  theory  {e.g.  Holton  (1979)),  and  includes 
Doppler  shifting  due  to  the  mean  environmental  flow. 

Given  the  area  bounded  by  the  gravity  wave,  the  individual  grid  box  tendencies  are 
computed  using  the  grid  volume  averaging  operator  deflned  by: 

G(x)  =  ffcXc  +  (1  -  o'c)®e,  (®) 

where  ac  is  the  fractional  area  within  a  grid  box  occupied  by  clouds.  In  our  case,  it  will  be 
the  total  area  occupied  by  the  active  population  divided  by  the  grid  box  area.  We  require 
that  the  sum  of  the  model  grid  tendencies  within  the  CIA  be  equivalent  to  the  average 
tendencies  computed  from  the  Arakawa-Schubert  parameterization.  Therefore,  if  we  let  Oij 
be  the  area  fraction  of  the  total  CIA  of  a  grid  box  (*,  j)  s.t  least  partially  located  within  the 
CIA,  any  average  quantity  must  conform  to  the  following: 

|=£o«G(eiy).  P) 

CIA 

The  model  grid  tendencies  are  computed  such  that  the  average  of  the  tendencies  in  the 
model  grid  box  which  contains  the  sub-ensemble,  and  the  tendencies  in  the  neighboring  grid 
boxes  which  are  within  the  CIA  are  exactly  equal  to  average  tendency  computed  over  the 
CIA.  Solving  the  system  of  equations  above  provides  the  grid  box  tendencies. 

Cloud  Propagation 

Each  sub-ensemble  may  advect  from  one  grid  box  to  another.  We  use  the  doud-base 
mass  flux  obtained  from  the  prognostic  mass  flux  equation  as  a  tracer  for  each  sub-ensemble. 
For  simplidty,  we  choose  the  advective  velocity  to  be  the  average  horizontal  wind  over  the 
depth  of  the  doud.  Eventually,  this  may  be  modified  to  indude  other  terms,  for  example, 
to  modd  the  effect  of  gust  fronts  on  storm  propagation,  or  storm  vortidty  or  helidty  to 
model,  right  and  left  moving  cells. 

SiiTTimary  of  Parameterization  Methodology 

We  now  summarize  the  methodology  of  the  parameterization: 

•  Using  a  cumulus  ensemble  model,  determine  the  sub-ensemble  properties  and  doud 
work  function. 

•  Integrate  the  prognostic  mass  flux  equation. 

•  Using  the  Simpson- Wiggert  method,  calculate  the  vertically-integrated  kinetic  energy 
of  each  sub-ensemble. 

•  Solve  for  the  quantity  ncr-the  total  sub-ensemble  area. 

•  Determine  the  PDF  and  select  an  active  doud  population. 

•  Determine  the  CIA  for  each  sub-ensemble. 

•  Determine  the  tendendes  valid  over  the  CIA  by  taking  the  divergence  of  the  doud 
fluxes. 

•  Determine  the  modd  grid  tendencies  such  that  they  average  to  the  CIA  tendency 
uRiTig  the  grid  volume  averaging  operator. 
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•  Advect  the  cloud  sub-ensembles  using  the  cloud-base  mass  flux  as  a  tracer. 

•  Repeat  the  procedure  on  the  next  time  step. 

1.2.4  Conclusion 

A  cumulus  paiameterization  suitable  for  a  wide  array  of  model  grid  scales  has  been 
described.  The  parameterization  is  fundamentally  different  than  presently  available  param- 
eterizations  because  of  its  grid  size  independence,  its  horizontally  non-local  nature,  and  use 
of  a  PDF  to  determine  the  active  population  lErom  the  sub-ensemble.  The  scheme  is  nu¬ 
merically  efladent  so  that  regional  scale  and  operational  numerical  models  can  potentially 
benefit  greatly  from  its  implementation.  Forecasts  obtained  from  this  parameterization  will 
represent  the  most  probable  doud  forcing  based  on  our  determination  of  the  PDF. 

We  have  just  begun  to  code  the  parameterization  and  as  of  yet  do  not  have  any  results. 
At  this  point,  the  code  is  being  tested  for  accuracy  and  effidency.  We  do  expect  a  /?  version 
of  the  parameterization  to  be  implemented  shortly. 

Also,  we  have  been  conducting  three-dimensional  cloud  ensemble  simulations  to  try  and 
better  parameterize  some  of  the  unknown  parameters  such  as  a  and  tj.  The  data  from 
these  simulations  has  not  yet  been  analyzed. 

Mr.  t^aflfin  is  expected  to  complete  the  preliminary  implementation  and  research  asso- 
dated  with  this  parameterization  within  the  next  year.  Progress  made  in  the  future  will  be 
described  in  his  dissertation  as  well  as  other  intermediate  reports  to  the  AFOSR/ AASERT 
contract  under  which  he  is  currently  funded. 
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